Atomic Bloch-Zener oscillations for sensitive force measurements in a cavity 
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Cold atoms in an optical lattice execute Bloch-Zener oscillations when they are accelerated. We 
have performed a theoretical investigation into the case when the optical lattice is the intra-cavity 
field of a driven Fabry-Perot resonator. When the atoms oscillate inside the resonator, we find that 
their back-action modulates the phase and intensity of the light transmitted through the cavity. We 
solve the coupled atom-light equations self-consistent ly and show that, remarkably, the Bloch period 
is unaffected by this back-action. The transmitted light provides a way to observe the oscillation 
continuously, allowing high precision measurements to be made with a small cloud of atoms. 
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When quantum particles in a potential lattice are sub- 
jected to a constant force F, they execute Bloch-Zener 
oscillations (BZOs) [1 with a frequency 

LJB=Fd/h, (1) 

where d is the period of the lattice. This behavior was 
first demonstrated [2 with electrons in semiconductor su- 
perlattices, where a DC electric field provided the force. 
However, rapid dephasing due to impurities [3] has pre- 
vented BZOs from becoming useful in solid state devices. 

Cold atoms in optical lattices have recently provided 
an alternative realization of BZOs H E H E El [10] 
in which long coherence times are possible. Initially, 
acceleration of the lattice induced the oscillations, but 
in subsequent experiments [71 [H [10], gravity provided 
the required force. In p!0], the BZO damping time was 
12s, allowing some 4000 cycles to be measured over 7s. 
With such long coherence times, cold atom BZOs be- 
come suitable for high precision measurements, for ex- 
ample to determine the fine structure constant a [4 , to 
measure gravity [10 , or to explore Casimir- Polder forces 
[TT] . In the experiments to date, it has been necessary to 
reconstruct the oscillations by making destructive mea- 
surements at a large number of different times, each mea- 
surement requiring a new cloud of atoms to be trapped, 
cooled, and loaded into the lattice. The process is la- 
borious and suffers from shot-to-shot variations in the 
initial cloud conditions. In this letter, we discuss how 
the measurement could be substantially improved by us- 
ing an optical cavity to enhance the interaction of the 
atoms with the light. We show how the light transmit- 
ted through the cavity can provide an in vivo observation 
of the BZOs and we assess the extent to which this per- 
turbs the motion of the atoms. Finally we consider the 
statistical sensitivity of the method and show that it can 
yield high precision in a single shot. 

Let us take the cavity to be a vertical Fabry-Perot res- 
onator illuminated by a laser, which makes a standing- 
wave light field inside, shown in Fig.jlja). Cold atoms in 
this lattice execute BZOs under the influence of gravity. 
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FIG. 1: Schematic of the proposed experiment, (a) A cloud 
of cold atoms is held in a standing-wave optical trap inside 
a vertical Fabry-Perot cavity, (b) The atoms execute Bloch- 
Zener oscillations, leading to a periodic modification of their 
wave function, (c) This modulates the intra-cavity power and 
hence the lattice depth s. (d) The power modulation is seen 
in the light transmitted by the cavity. 



as illustrated in Fig.jTJb). The optical dipole interaction 
between one cavity photon and an atom placed at an 
antinode of the field is given by %o = f^\/^c/{^oV), 
where cOc and V are the frequency and volume of the rel- 
evant cavity mode and /i is the atomic transition dipole 
moment. The effect of one atom on the cavity field is 
characterised by the cooperativity C = ^q/(2a>:7), where 
27 is the atomic spontaneous emission rate in free space 
and 2/^ is the cavity energy damping rate. When C > 1, 
the cavity field is strongly perturbed by the atom, as 
illustrated in Fig.jljc). Thus the light transmitted or re- 
flected by such a cavity can detect the presence of a single 
atom [12^, and can be sensitive to the motion of atoms 
trapped within the cavity [13], as in Fig.[l];d). Although 
we shall not discuss Bose-Einstein condensates (BECs) in 
this paper, we note in passing that several experiments 
have succeeded in placing BECs inside optical cavities 
[14. 

Consider a cavity mode, whose frequency in the ab- 
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sence of any atoms is c^c, pumped by an external laser of 
the same frequency. A cloud of N atoms, each of mass 
m, is placed inside the cavity and is sufficiently dilute 
that the atoms do not interact directly with each other. 
The coupled atom-cavity hamiltonian becomes [15] 
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where the operator d{t) annihilates a photon in the cav- 
ity mode and the operator il){z^t) annihilates an atom 
at the point z. The first term gives the kinetic energy 
of the atom. The second describes the quadratic Stark 
interaction in the rotating- wave approximation. Here, 
A = — is the detuning between the cavity mode 
and the atomic transition frequency a;^, and kc = oOc/c. 
This is an approximate form, that is valid when A 7 
and T]'^ <C (hzA/go)'^. Under these conditions the atom 
has negligible population in the excited state. The third 
term accounts for the external force, and the last term 
describes the coherent excitation of the cavity by the ex- 
ternal laser. For a cavity with equal mirror reflectivities, 
the pumping rate is 7^ = where / is the rate of in- 

cident photons matching the cavity mode. The cavity 
field is a driven and damped quantum harmonic oscil- 
lator for which it is known that exact solutions of the 
Fokker-Planck equation are coherent states [16 . Fur- 
thermore, since there is negligible spontaneous emission 
by the atoms, the cavity field remains in a coherent state 
in the presence of atoms. Taking the expectation value 
of the Heisenberg equations of motion for d and lip in the 
coherent state \a) yields the equations of motion pT] : 
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where "^{z.t) = (^Ar-i|^(^, ^)|^Ar)/vA^ is the wave 
function occupied by all N atoms and I^at) is the cor- 
responding state vector in Fock space. We have added a 
damping term proportional to k, in Eq. ([3| to account for 
leakage of light through the mirrors [18^ . These equations 
neglect quantum fluctuations of the light field which can 
heat the atoms: we return to this effect later. Equa- 
tions ([3| and Q must be solved self-consistently: the 
coupling 
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changes a, which changes the depth of the lattice. This 
alters the atomic wave function and therefore changes 
etc. In static equilibrium, a = and then 

(6) 



a = 
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Even if the atoms are in motion, ^ remains a very 
good approximation since is generally much greater 
than the highest frequency in the atom dynamics, so that 
the field 'instantaneously' adapts to the atomic distribu- 
tion. Frequencies that feature in the atomic motion are 
the BZO frequency, the band splitting, and the harmonic 
frequency cJho = 2^o|<^| V^r/ (^^) at the bottom of each 
potential well, E'r = li?k'^/{2m) being the atomic recoil 
energy. For the experiments we consider here cjb and cjho 
are much smaller than so we assume in our analytic 
calculations, though not in our numerical simulations, 
that Eq. ^ holds. 

Let us recall the standard theory of BZOs without 
a cavity [El |20l ED EH |23l El |25j. The Schrodinger 
Eq. Q can be written as ifi^ = (^0 + Fz)^ ^ where 
^0 = /{2m) + 5cos^ kcZ. The eigenfunctions Xg,s,n(^) 
of Hq are Mathieu functions (Bloch waves) that in gen- 
eral depend on position lattice depth 5, band index n 
and quasimomentum (7, restricted to the first Brillouin 
zone —Ti/d < q < n/d [26 . We assume in our analytic 
calculations that the atoms are in the lowest band, whose 
energy is £^g,s, and dispense with the band index so that 
HoXq,s = Eq^sXq,s' The Bloch theorem allows us to write 
Xq,s{z) = Uq^s{z)exp[iqz], where Uq^s{z) = Uq^s{z ^ d) 
obeys 
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Uq^s{z)^SCOs'^{kcZ)Uq^s{z) = Eq^sUq^s{z). (7) 



To tackle the full hamiltonian Hq-\-Fz we make the gauge 
transformation "^{z^t) = exp[— iFtz/fi]^(z, t), yielding 

the Schrodinger equation ih^ = where H = {p — 

FtY/{2m) + 5cos^/Cc^- Comparing this with ^ we 
see that the effect of the force is to evolve the quasi- 
momentum according to Bloch's acceleration theorem pLj 



q qit) = qo- Ft/h, 



(8) 



where go is the quasimomentum at t = 0. When q{t) 
reaches the edge of the Brillouin zone at —ir/d it is 
mapped to the identical point q = +7r/(i, giving rise to os- 
cillatory behaviour - the BZO. The corresponding Bloch 
wave has the approximate form [27 (setting qo = 0) 



^{z,t) ^ /7,(,),,(z)exp[-i/n / dt' E^ 



(9) 



l + iAr^2(t)/(/^A)* 



within the adiabatic approximation that the rate of 
change U/U is too small to excite higher bands. Here 
Uq{t),s{^) is the instantaneous solution of Eq. Q. Dur- 
ing a BZO the spatial distribution Uq(^t),s{^) oscillates 
with a breathing motion, as shown schematically in Fig. 

1(b). 

Consider now the effect of the BZOs on the field inside 
the cavity. The coupling g (Eq. ^) depends on |^(z,t)p 
which equals \^{z,t)\'^. Its breathing motion changes g, 
which in turn modulates the cavity field through Eq. (pi) . 
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Inserting (|6|) into (|4|), and replacing r] by Vi^I, we ob- 
tain the Schrodinger equation for atoms in a periodic 
potential ih^f = {p^ /2m + s(t) co^^ (k^z) + Fz)^, with 
the time-dependent potential depth 



s(t) = hi 



Ak 



{NgHt)/{AK)y 



(10) 



When NgQ/{hzA) <C 1, 5 is approximately constant in 
time and a becomes {r]/hz)[l — iNg{t)'^ /{Ahz)], i.e. the 
light exhibits a small phase modulation and has negligi- 
ble intensity variation. In this case, the atoms oscillate 
in a lattice that is essentially static. With stronger cou- 
pling, where NgQ/{HiA) :^ 1, s{t) is changed significantly 
during an oscillation. The fundamental period is never- 
theless unchanged and is still given by Eq. 0. Physi- 
cally, this is because BZOs arise from an interference of 
waves in a lattice akin to Bragg scattering and lattice 
depth plays no role in determining the phase-matching 
condition. Rather, this is determined by the symmetry 
of the hamiltonian which is precisely maintained at all 
times. To examine the effect of lattice depth modula- 
tion in more detail, let us begin with the adiabatic case 
where the frequency spectrum of s{t) remains largely at 
low frequencies unable to excite higher Bloch bands. An 
example is shown in Fig. [2|a). In that case, the low- 
est band energy £^q,s(t) is still determined by the instan- 
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and qo remains a constant of the motion generated by 
Ho{t) despite the time- dependent potential ([19] reaches 
a similar conclusion for electrons in an ac field). Con- 
sequently, the Bloch wave in the presence of an external 
force can still be calculated using Eq. ([9|, provided we 
use the instantaneous values of s{t) and s{t'). It only re- 
mains to find the self-consistent solution for s{t) by solv- 
ing Eq. (10) at each instant of time. Here s{t) appears 
explicitly on the left and also implicitly on the right as 
a parameter determining the Bloch wave function that is 
required to calculate the coupling \/Ng(t) using Eq. ([s]). 
We conclude that despite the lattice depth modulation, 
the Bloch acceleration theorem (|8| still holds for atoms 
in the lowest band and therefore the fundamental oscilla- 
tion frequency remains identical to the cjb of an atom in 
a static lattice. Furthermore, because cjb is the same for 
all bands, there is no frequency shift even when higher 
bands are excited, as we have verified numerically for a 
wide range of NgQ/{Ahz). 

Fig. |2|a) compares the lattice depths obtained by nu- 
merical solution of Eqs. ([3| and Q (lines) and by the adi- 
abatic approximation of Eqs. (|9f and ([lO| (dots). The 
two are in good agreement. Figure[2|b), shows a case 
of stronger coupling, where one can see fast oscillations 
superimposed on the main motion. The Fourier trans- 
form of this reveals two effects, illustrated in Fig.[2|c). 
(i) There are higher harmonics of cub because the os- 
cillations at the Bloch periodicity are not exactly sinu- 
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FIG. 2: Calculated evolution of lattice depth s{t) normal- 
ized to the recoil energy for ^^Rb atoms undergoing 840 Hz 
Bloch-Zener oscillations in a 780 nm lattice with s ~ 3^^. 
Lines: numerical solution to Eqs. ([3| and p]). Dots: self- 
consistent adiabatic approximation of Eqs. (joT and (10). (a) 
^9o/{i^A) = 0.4. (b) Close-up of lattice depth oscillations 
with stronger coupling, NgQ/{K,A) = 1. Non-adiabatic ef- 
fects are seen in the line, (c) Fourier transform 5(0;) of result 
in (b), showing harmonics of fundamental frequency ujb- In- 
set: Close-up of s{oo) at higher frequencies. Harmonics of oob 
appear as sharp vertical lines. In addition, one sees much 
weaker, broad, Fourier components due to band excitation, 
corresponding to the rapid oscillations in (b). The actual in- 
dividual values of the parameters used (or assumed in the case 
where only ratios enter) in the calculations were (see text): 
N = 5 X 10^, ^0 = 27r X 2.8 MHz, n = 27t x 1.0 MHz. In 
(a) A = 27r X 1.0 THz, ry = 27r x 39 MHz; in (b) and (c) 
A = 27r X 0.39 THz, ry = 27r x 28 MHz. 



soidal. In the adiabatic solution, the first four are ac- 
curately reproduced and the higher harmonics are very 
small, (ii) The exact solution of Eqs. (|3| and (|4| predicts 
non-adiabatic components at higher frequencies, shown 
inset in Fig.[2|c). These are predominantly harmonics 
of cob (the sharp lines), but in addition, there are other 
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frequency components that can be seen as broad lines. 
These are only found in the full numerical solution and 
are due to a small amount of excitation to higher Bloch 
bands. These non-adiabatic effects become stronger as 
the parameter NgQ/^n/S) is increased. 

The atoms are driven not only by the mean-field po- 
tential s{t) cos^(/Cc2:), but also by random forces due to i) 
spontaneous emission, and ii) fluctuations in the photon 
number, that are associated with the decay rates 7 and 
respectively. In particular, in the strong coupling regime 
photon number fluctuations can signiflcantly heat atoms 
inside an optical cavity [13] . We quantify the heating ef- 
fect via the increase in the width dp of the atom's momen- 
tum distribution according to da^/dt = 2D. The diffu- 
sion constant D = dt'[{F^^(t)F^^(t + t')) - {F^^)^] 
[28] involves two-time correlations of the dipole force 
Fdip, and hence of the intracavity electric fleld. It has 
been calculated for atoms in cavities in [29] and has two 
terms. The flrst occurs in any standing-wave light fleld 
and at low saturation is given by I^sw = ^^^^/2rsp [28 , 
where = 2^\a\^gQ//S? is the spontaneous emission 
rate at an antinode. The second term, speciflc to cav- 
ities, is I^cav = 2I)sw C'sin^(2/cc^). This diffusion lim- 
its the coherent measurement time, which we take to be 
the time r when the momentum distribution has a width 
equal to one half of the flrst Brillouin zone, i.e. <jp = hkc 
Then r = rsp/(l + where we have replaced s\v?{2kcz) 
by 1/2 - a good approximation in the ground band for a 
lattice of depth s = 3E^. 

The BZOs can be observed by detecting the photon 
current \a{t)\'^hz transmitted through the cavity which 
is directly proportional to the depth s{t) of the lattice 
{s = %ol<^P/^5 Eq. Q), whose evolution is shown 
in Fig. [2] For an estimate of the measurement preci- 
sion, let us write the detection rate as R[l + e cos(cjt)]. 
After measuring this for a time r with detectors hav- 
ing an efliciency ^, the shot noise gives an uncertainty 
in the oscillation frequency of ^ 2'kt~'^^'^ / {cy/^R)^ 
in which comes from the linewidth due to the fl- 
nite duration of the measurement, and r~^/^ comes from 
the shot noise in this bandwidth. This simple estimate 
is close to the Cramer- Rao lower bound [30 , the limit 
given by the information content of the signal. For 
small e, R ~ lap/^:, then the frequency uncertainty can 

be written as cr^ ^ ^^i^i^)^ih + ^^^^^ 
to bring out the implicit dependence on the number of 
atoms N contained in this result, we deflne the parame- 
ter X = NgQ/{hzA). Referring to Eq. (10), if the number 
of atoms is increased then proportional increases in the 
laser detuning and intensity maintain constant values of 
5, X and e, while the measurement time r and the intra- 
cavity power both increase by the factor N. With this 
scaling, the measurement time, detuning and frequency 



uncertainty are given by 
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(12) 
(13) 



The uncertainty cr^ therefore decreases rapidly with a 
dramatic dependence which ultimately derives from 
the continuous observation of the oscillations (via the 
enhanced measurement time r). 

Small s reduces a^^ but the lattice must support the 
atoms against gravity. We flnd that s = 3Eji is a rea- 
sonable compromise. In the example of Fig.[2|a) we have 
chosen x = 0.4, which gives e = 1.3% in this lattice. 
Taking a reasonable number of atoms, let us say 5 x 10^, 
a readily achieved cooper at ivity of C = 1.3, and a pho- 
ton detector with 60% efficiency, brings (Tuj/^b to Ippm. 
From the deflnitions of x and C this requires a detun- 
ing of A = 27r X 1 THz and from Eq. (11) a measure- 
ment time of only r = 1 s. These numbers flx the ratio 
g'^/tz = 27r X 7.8 MHz. If we choose 1^ = 27t x 1 MHz, 
then using s = %ol<^P/^5 this means there are on aver- 
age > 1400 photons in the cavity. 

We have not included direct atom-atom interactions 
in the model discussed here. These can lead to 
quasimomentum-changing transitions (non- vertical tran- 
sitions in the language of [31]) which dephase the BZOs. 
However, as summarized in the opening paragraphs of 
this paper, long-lived BZOs have already been success- 
fully demonstrated in cold gases containing many atoms, 
and so it is a question of degree, i.e. at what atom den- 
sity and interaction strengths do the interactions be- 
come important? An experiment investigating the con- 
trol of interaction-induced dephasing of BZOs in a Bose- 
Einstein condensate has recently been reported [32] . Us- 
ing a Feshbach resonance they were able to increase their 
dephasing time from a few to more than 20 thousand 
BZO periods. From the details of their measurements we 
estimate that, for our example given immediately above 
involving 5 x 10^ atoms, the dephasing due to collisions 
is negligible for reasonable cavity geometries. The effect 
becomes signiflcant on increasing the number of atoms 
to several million, but can be suppressed by tuning to a 
Feshbach resonance [32l . Large detuning and laser power 
impose a practical limit on the useful atom number at 
about this level anyway. 

For atoms being continuously measured, an impor- 
tant source of dephasing is quantum measurement back- 
action. This effect is included in the estimate above in 
a quasi-classical way through the diffusive heating of the 
atoms by fluctuations of the cavity light fleld. The cavity 
fleld suffers fluctuations because it is dissipatively cou- 
pled to the outside world and it is precisely the light 
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escaping from the cavity (at rate k,) that contains the in- 
formation about the state of the atoms. In other words, 
it is the cavity decay that is doing the measuring. For 
our parameters, r < ^Tgp because C > 1, and therefore 
it is the fluctuations due to cavity decay that hmit the 
measurement time, rather than the spontaneous emis- 
sion. We plan to perform a more microscopic study of 
the measurement back-action in the future. 

In conclusion, we predict that the force on a small cold 
atom cloud can be measured very accurately by a new 
method based on BZO oscillations in an optical cavity. 
The BZO oscillations are measured continuously by mon- 
itoring the light that leaks out of the cavity. This enables 
a relatively fast and high precision measurement of the 
oscillation frequency, the error being given by Eq. (13). 



Our treatment of the problem is based upon solving the 
coupled equations of motion for the atoms and light. This 
gives a detailed picture of the dynamics, including the 
adiabatic and non-adiabatic aspects. 

We note that since the submission of this paper a re- 
lated proposal on monitoring of Bloch oscillations using 
an optical cavity has appeared [33 . 
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